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Abstract 

There are recent interests with CsI(Tl) scintillating crystals for Dark Matter 
experiments. The key merit is the capability to differentiate nuclear recoil (nr) 
signatures from the background /3/7-events due to ambient radioactivity on the 
basis of their different pulse shapes. One of the major experimental challenges is to 
perform such pulse shape analysis in the statistics-limited domain where the light 
output is close to the detection threshold. Using data derived from measurements 
with low energy 7's and nuclear recoils due to neutron elastic scatterings, it was 
verified that the pulse shapes between /3/7-events are different. Several methods of 
pulse shape discrimination are studied, and their relative merits are compared. Pull 
digitization of the pulse shapes is crucial to achieve good discrimination. Advanced 
software techniques with mean time, neural network and likelihood ratios give rise 
to satisfactory performance, and are superior to the conventional Double Charge 
method commonly applied at higher energies. Pulse shape discrimination becomes 
effective starting at a light yield of about 20 photo-electrons. This corresponds to a 
detection threshold of about 5 keV electron-equivalence energy, or 40—50 keV recoil 
kinetic energy, in realistic experiments. 
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1 Introduction 



The detection of Dark Matter and the studies of their properties are of fundamental 
importance in particle physics and cosmology. The Weakly Interacting Massive Particles 
(WIMPs) are good candidates for "Cold Dark Matter", and their experimental searches 
have gathered a lot of interests in recent years. The most promising avenue is to detect 
the nuclear recoil signatures due to elastic scatterings of WIMPs on the target isotopes. 
The typical energy depositions are only of the order of 10 keV, imposing big experimental 
challenges in terms of the detection of weak signals as well as background control at low 
energy close to detection threshold. A wide spectrum of experimental techniques is being 
pursued P|. There is still much room for new detector concept to push the sensitivities 
further. It would be of great interest if the sensitivities of WIMP searches can probe the 
level predicted by the various Super-Symmetry models. 

There are potential merits of using CsI(Tl) scintillating crystals [5^ for WIMP search 
and other low-energy low-background experiments jSl 13 • An experiment with 200 kg of 
CsI(Tl) crystal scintillators to study low energy neutrino interactions at the Kuo-Sheng 
power reactor is being pursued |11E|, while the adaptation of the crystal for Dark Matter 
searches is the focus of several R&D projects [HI 13 IS] and an on-going experiment P]. 

The high-A content of the Csl enhances the sensitivities for the spin-independent 
interactions (which depends on the neutron number squared) between the WIMPs and the 
target, relative to most other candidate target isotopes. The high-Z composition allows a 
compact design and provides large suppression of background due to ambient radioactivity 
if a three dimensional fiducial volume definition can be realized. Both ^^^Cs and ^^^I are 
100% in their respective isotopic abundance. Being close in their mass numbers, the 
response to nuclear recoil from the interactions with WIMPs would be similar, allowing 
simpler interpretations of the experimental signatures. 

As a detector, the crystal has large light yield, low energy threshold and with pulse 
shape discrimination (PSD) characteristics for differentiating /3/7 background from the 
nuclear recoil events Scintillating Nal(Tl) crystals with the order of 100 kg target 

mass have been deployed for Dark Matter experiments JU]; but it has been shown that 
CsI(Tl) provides superior PSD capabilities to Nal(Tl) Unlike Nal(Tl), CsI(Tl) is 

only slightly hygroscopic such that it can be machined easily and does not require hermetic 
seal (that is, passive materials) in a large detector system. In addition, large (40 tons) 
electromagnetic calorimeter systems ^1] have been constructed and made operational in 
high energy physics experiments, making this technology affordable and realistic to scale 
up. Considering all the associated costs, the price of CsI(Tl) is in fact less than that for 
Nal(Tl). In order to produce positive and definite evidence of the WIMPs, an accurate 
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measurement of the annual modulation (where the maximal effects are only 7%) would 
be necessary such that the availability of large target mass is a very desirable feature. 

One of the key issues to realize a Dark Matter search experiment with CsI(Tl) crystal 
scintillator is the studies of the experimental signatures of nuclear recoils due to WIMP- 
nuclei elastic scatterings. Nuclear recoils produce high charge density (DA/dx) such that 
the scintillating light yield is "quenched" and the timing profile of pulse shape is different 
relative to the same energy deposition by minimum ionizing particles ^21 • These WIMP- 
induced signatures are the same as the nuclear recoil events produced by elastic scattering 
of neutrons on nuclei, and hence can be studied in the laboratory. 

2 Pulse Shape Discrimination 

It has been well-studied [21 E] that the light emission profiles of scintillating CsI(Tl) crys- 
tals exhibit different shape for 7-rays and electrons (that is, minimum ionizing particles), 
as compared to that for a-particles and nuclear recoils at the high energy (>100 keV) 
regime . Heavily ionizing events due to a-particles and nuclear recoils have faster de- 
cays than those from e/7's — opposite to the response in liquid scintillator This 
characteristic property makes particle identification possible with this scintillator [TH] . 

Matured pulse shape discrimination (PSD) techniques have been devised at high en- 
ergies where the photo-electrons are abundant. The experimental challenge for adapting 
the PSD idea to Dark Matter experiments is that one must now work in the regime where 
the number of photo-electrons (Npe) is small such that the statistical fluctuations may 
wash out the differences. 

In the following sub-sections, we verify that the pulse shapes of 7 and nuclear recoil 
events with CsI(Tl) crystals are different even at the Dark Matter relevant low energy 
regime (<100 keV). The theme and focus of this article is to investigate and compare 
the various software techniques which can perform event-by-event PSD at this low light 
output domain. Detailed response characteristics of CsI(Tl) crystals at the low energy 
regime have already been studied in previous work [ni[7[|HlE]- 

2.1 Measurements 

A CsI(Tl) crystal of dimensions 5 cm x 5 cm x 5 cm and mass 560 g was used to provide 
data for these investigations. The light emissions were read out by a 29 mm diameter 
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photo-multiplier tube (PMT)^ with standard bi-alkali photo-cathode. The conversion 
factor between energy deposition and light output is 4 photo-electron per keV of electron- 
equivalence (e.e.) energy. This was obtained by calibration measurements with an LED 
pulser operated at the Npe ~ 1 intensity. The events were digitized by a 20 MHz (that 
is, 50 ns for one time-bin) Flash Analog-to-Digital Convertor (FADC) ^3] with 8-bit 
resolution, such that the pulse shape can be denoted by an amplitude time-sequence Aj. 

The low energy 7-data were taken with a standard radioactive ^^^Ba source, which 
provides several 7-lines up to 356 keV. The series of associated 7-peaks provide good 
energy calibration at low energy. The low energy (<30 keV) events crucial for this study 
are due to Compton scatterings of the higher energy 7's such that they originate from the 
bulk of the crystal. These merits justify the choice of ^^^Ba source over other low energy 
sources like ^^Fe (5.9 keV) and ^°^Cd (22.1 keV) where the attenuation of the 7's due to 
the crystal wrapping materials is severe, and the events only originate at the surface of the 
crystal. Nuclear recoil data, on the other hand, were taken from the neutron facility at the 
13 MV Tandem accelerator at the China Institute of Atomic Energy at Beijing. The data 
consisted of Time-of-Flight (ToF) measurements which helped to distinguish nuclear recoil 
from the other background events. The results of the quenching factor measurements were 
already published jHj. For completeness, a-events with ^^^Am (5.49 MeV) source were also 
recorded. 

The nuclear recoil pulses recorded in a neutron beam environment were contaminated 
by an intense accidental 7-background. The average pulse shapes for both nuclear recoil 
and 7-background (as identified by the ToF cut) events derived from the neutron beam 
measurements at a nuclear recoil e.e. energy of 4.8 keV are depicted in Figure H The long 
tails indicate there is a substantial contribution from time-uncorrelated 7-background. 
Upon taking averages from a large sample and subtracting the 7-background, such data 
are sufficient to provide a good quenching factor measurement as well as the average 
"background-free" nuclear recoil pulse shapes, as depicted in Figures Et&b. 

2.2 Average Light Emission Pulse Profiles 

From the FADC measurements with 7-events and the subtraction procedures for recoil 
events discussed in the previous section, the average pulse shapes for both categories are 
depicted in Figures EK&b. 

In Figure ISt', the spread among the average 7-pulses between 5 keV and 40 keV are 
denoted by dotted lines, as compared to the solid line for the recoil pulse shape at 4.8 keV 
^CRllO, Hamamatsu Photonics, China 
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e.e. energy (43 keV recoil energy). Similarly, the spread among recoil pulse shapes between 
4 keV and 11 keV e.e. energy (about 40 keV to 110 keV recoil energy) are defined by the 
dotted lines in Figure |2b; whereas the solid line denotes 7-profiles at 20 keV. It can be 
concluded that the energy dependence of the pulse shapes within the recoil and 7-data 
samples at this low (<100 keV) energy range is small compared to the differences in the 
pulse shapes between the recoil and 7-events. These results are consistent with those from 
Refs. 13 13 El • Event identification is in principle possible be achieved at low energy with 
CsI(Tl) crystals from the pulse shape information. 

The pulse shapes(A) as a function of time(t) displayed in Figures Efe^&b can be fitted 
to an analytical form 

t 1 t r t 

A = Constant * [ 1 — exp( ) ] * [ — exp( ) H exp( ) ] , (1) 

To Ti T~i T2 T2 

where tq is the rise time, {ti,T2) denote the fall times, and r is the ratio between the 

slow and fast decay components. As illustrations of the typical ranges, the best-fit values 

for recoil events at 4.8 keV e.e. energy and 7-events at 20 keV are tabulated in Table ^ 

Those for undoped Csl crystal are also shown for comparison. The values of tq in CsI(Tl) 

are dominated by the electronics shaping rise time of 250 ns for > /is pulses jli]. The 

intrinsic rise times of the CsI(Tl) scintillator are expected to ~125 ns and ~20 ns for 7- 

and recoil-events, respectively The difference in the decay time constants between the 

recoil and 7-events is the basis of pulse shape discrimination. 

2.3 Data Samples 

The direct measurements of Figure ^ indicate that the average recoil pulse shape can 
be derived by statistical subtraction. However, at the event-by-event level, the time- 
profile for photo-electron emissions of the neutron beam data set is complicated by an 
uncontrolled and sizable background contribution. Dark Matter searches, on the other 
hand, are low-count-rate experiments such that nuclear recoils due to WIMP interactions 
will not be contaminated by accidentals. Therefore the neutron beam data do not provide 
a realistic sample for the studies of detector response and realistic signals in WIMP 
searches at the event-by-event level. 

As remedies, the single-event nuclear recoil pulse shape was generated by simulations, 
where the input included the measured pulse profile at 4.8 keV e.e. energy from FigureEK 
and the parametrization of Tabled The simulated events were convolutions of (a) a total 
of Npe single photo-electron pulses whose timing fiuctuations were generated according 
to the average recoil pulse shape, and (b) the single photo-electron response of the PMT 
and readout system for each of these pulses, provided by the LED pulser measurements. 
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A "self-trigger" criterion was imposed to mimic the realistic situation — that is, the time- 
zero of the events was defined by the first instants where the pulse was above a specified 
threshold. As illustrations, typical events at Npe ~ 20 from the measured 7 and simulated 
nuclear recoil data are displayed in Figures OK&b, respectively. Both categories of events 
are similar by visual inspection, demonstrating that (a) the simulation algorithms are 
valid, and (b) advanced pattern recognition techniques would be necessary to achieve 
event identification. 

Applying the same algorithm on the 7- reference profile at 20 keV shown in Figure Eb, 
simulated 7-events were also generated. The comparisons between the distributions of the 
PSD figures-of-merit among the simulated and measured 7-events at Npe ~ 60 are depicted 
in Figures E^&b. The agreement is excellent, further demonstrating the ability of the pulse 
shape simulation algorithms in reproducing correctly the fluctuations among individual 
events. The figures-of-merit will be defined and discussed in details in Section 12.41 

The excellent consistencies in Figures E^&b between simulations and data for the 
Npe ~ 60 7-events justifies the use of the simulated nuclear recoil and 7 data set (denoted 
by Dnr and D^, respectively) for the PSD studies discussed in subsequent sections. There 
are better uniformity and systematic control among the different simulated data set such 
that the residual systematic effects will be minimized and canceled out when comparisons 
are made. The input parameters to the simulation procedures can be varied to study 
features like sensitivities and robustness. 

2.4 Classical Pulse Shape Discrimination Method 
2.4.1 Double Charge Method 

A well-established way to achieve PSD at high light yield is the "double charge method" [T5] 
This involves the comparison of the "total charge" (Qt) and the "partial charge" (Qp), 
which are the total and partial integration of the pulse, respectively. This is the standard 
approach with Analog Digital Convertor (ADC) based data acquisition systems where the 
complete pulse shape information is not available. Typically, the partial charge measure- 
ment is done by delaying the PMT pulses via cabling and both the prompt and delayed 
signals are read out by the ABC sampled with the same gate. 

Displayed in Figure is the comparison of 7 and a events at the MeV energy range 
from data with ambient radioactivity and ^"^^Am a-source, respectively. The ranges were 
chosen such that Qt and Qp involve integration over 4 //s after trigger and after a delay 
of 0.5 fis, respectively. A j/a separation of >99% efficiency down to about 200 keV e.e. 
light output can be achieved. It has been shown that PSD can be achieved even in high 
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energy events where the FADC measurements are saturated PF. However, as indicated 
in Figure El one would come into difficulties to perform PSD with this simple algorithm 
for events at light yield below 100 keV e.e. energy. 

2.5 Pulse Shape Discrimination Methods at Full Digitization 

With the advent and popular usage of FADCs, complete pulse shape digitization becomes 
realistic. Three different pattern recognition techniques were investigated, all of which 
rely on the full digitization of the PMT signals. 

2.5.1 Mean Time Method 

The measurement of the average time for individual events by the mean time (MT) method 
has been used for PSD studies p]. The mean time is defined as 

E (Ai ti) 

(t) ^ . (2) 

i 

where Ai is the FADC-amplitude at time-bin tj. 

The typical (t) distributions at Npg ~ 20 for Dm- and D^ are displayed in Figure IHt', at 
an integration of 5 fis after the time-zero set by the trigger. It can be seen that satisfactory 
separation can be achieved under such conditions. 

2.5.2 Neural Network Methods 

The neural network (NN) methods J7j are now frequently adopted for analysis in high 
energy physics experiments. It has been applied to event-by-event pulse shape analysis 
for background identification in double beta decay searches jTHj. The pedestal-subtracted 
FADC data within 5 fis after trigger corresponds to the input nodes of the neural network. 
That is, the network has Ni=100 input nodes denoted by X(xi) with the integrated sum 
normalized to unity: 

Ni 

E^i = l- (3) 

i=l 

Negative values were reset to zero. In addition, the number of hidden nodes was selected 
to be Nh=25. It has been checked that the results are independent of this choice, so long 
as Nh > 20. 
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Adopting the Neural Network JETNET 3.0 package J7], a function F(X) is defined 
such that 

Nh N; 

F(X) = G(^ Uj G(^ wjk Xk + ^^j) + 0o) (4) 
j=i k=i 

where (uj, wjk) and (^j,0o) are the "weight" and "offset" coefficients, respectively, to 
be derived from the training samples, and the function G(y) is the non-linear neuron 
activation function 

G(y) = ^ [ 1 + tanh(y) ] = , (5) 

which is the functional form characterizing a 3-layer neural network consisting of the 
input, hidden and output layers. 

A total of Nt=4000 events from both the D^r and data sets are used as train- 
ing samples, corresponding to T(X)=1 and 0, respectively. The optimal coefficients are 
obtained by minimizing the error function 

Nt 

E = E [F(X)-T(X)]2 . (6) 

i=l 

Once the coefficients are derived, the procedures are applied to independent recoil and 
7 data set. The F(X)- values for recoil events would be larger than those for 7-events at the 
same light yield. The comparisons of F(X) distributions for the simulated and measured 
7-data at Npe ~ 60 are shown in Figure |^, while those for Dnr and at Npe ~ 20 are 
displayed in Figure Eb- It can be seen that the simulations agree well with data, while 
there is good separation between the recoil and 7 samples. 

2.5.3 Likelihood Ratio Methods 

Motivated by the commonly-used of likelihood ratio test PIE] for the goodness-of-fit, a 
likelihood ratio (LR) method was devised to perform the tasks of pulse shape analysis. 
Similar methods are successfully applied in high energy physics data analysis in comparing 
likelihoods and assigning probabilities among the different hypotheses for events where 
many output parameters are measured. The reference profiles for neutrons and 7's from 
Figure |21 are required as the input. This is different from the previous two techniques 
where prior knowledge of the reference profiles is not necessary. 

The areas of the reference pulses are normalized to unity, and the profiles are denoted 
by arrays R(ri) and r(7i) for the nuclear recoil and 7 reference shapes, respectively. Two 
likelihood functions, Lj. and L^, are defined for each event: 

N; N; 

Lr = n ; L7 = n 7r , (7) 

i=l i=l 
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where X(xi) with dimension N; = 100 are the measured pulse shape information for the 
events to be analyzed, as defined in Section 12.5.21 The likelihood functions quantify how 
probable the measured pulse shapes do originate from the the reference profiles. The 
likelihood ratio LR defined by: 

^ l:tl; • 

will test which hypothesis is more likely. The LR-values for recoil events would be larger 
than those for 7-events at the same light yield. 

The comparisons of LR between simulated and measured 7-data at Npe ~ 60 are 
shown in Figure |3Jd. while the typical LR distributions for the Dm- and data set at 
Npe ~ 20 are depicted in FigurelHfc. Similar to the neural network methods, there are good 
agreement between simulations and data, while there is satisfactory separation between 
the recoil and 7 samples. 



2.6 Comparisons 

The excellent agreement depicted in Figures E^&b between measured and simulated 7- 
data set justifies that valid comparisons can be made between the simulated data set 
and Djir- To quantify, two figures of merits are defined: (a) ego: the survival efficiencies of 
Dnr at selections which ensure that 90% of the events are suppressed; and (b) /90: the 
probabilities where the events would be mis-identified as recoil signals at cuts where 
90% of would survivc. Both ego and 190 are energy dependent, and would approach 1 
and 0, respectively, at the high light yield (large Npe) limits. 

The variations of ego and /90 as a function of Npe for the three different methods (MT, 
NN, LR) are depicted in Figures [7^&b, respectively. The photo-electron number Npe was 
adopted as the unit to characterize the light yield. In this way, the results can be directly 
applicable to other configurations using CsI(Tl) as the detector medium. Dotted lines 
in Figures [Tti&b corresponds to the survival probabilities of and Dnr, respectively. 
The results indicate that the generic features that all the three methods: (a) can achieve 
PSD with satisfactory efficiencies (>50% 7-background rejection) at Npe > 20; (b) can 
identify >90% of the D^ background while keeping the efficiencies for Dm- to be >90% at 
Npe > 80; and (c) give similar performance at the large light yield limit (Npe = 120), which 
approaches the expected values of 1 and for ego and Z90, respectively. Comparing the 
relative merits among the three algorithms, the NN technique gives better performance 
in the low energy range. The relative merits between the LR and MT methods are 
similar at low statistics (Npe < 20), while LR tends to perform better at the intermediate 
range (40 < Npe < 80). Accordingly, the neural network method would be the preferred 
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technique in pulse shape analysis at regime where the statistics is marginal. 

Tests have been performed on simulated events with different single photo-electron 
response functions. Results consistent with the performance parameters shown in Fig- 
ures Eti'&b were obtained. This shows that the results are robust and insensitive to the 
details of the simulation algorithms so long as the same reference profiles in Figures |2ti^&b 
are used for the photo-electron timing distributions. This indicates that the results would 
also be valid in measurements of CsI(Tl) crystals where the PMT response and electronics 
settings (such as shaping times) would be different. 

3 Summary and Conclusions 

This article reports the measurement of light emission pulse profiles for nuclear recoil 
and 7-events in CsI(Tl) crystal scintillator at the energy range relevant for Dark Matter 
searches (< 100 keV). The energy dependence of the profiles within the recoil and 7- 
samples is small compared to the differences between them. Event identification is feasible 
for CsI(Tl) based on pulse shape discrimination. 

Various software techniques to achieve pulse shape discrimination in this "near thresh- 
old" regime were studied. The performance of the three methods based on complete pulse 
shape information (mean time, neural network and likelihood ratio) is superior to the ma- 
tured and conventional double charge method well-demonstrated when photo-electrons are 
abundant. Full digitization is crucial for achieving PSD at the marginal statistics domain. 
Among the three methods studied, the neural network technique provides the best per- 
formance in regime where the statistics is marginal, while the other two methods are still 
satisfactory. The algorithms are robust and insensitive to the measurement parameters 
like PMT response or electronic shaping times. 

The results from this study are relevant to the potential capabilities and practical 
design of Dark Matter experiments based on the CsI(Tl) crystal. Satisfactory (bigger then 
50%) separation between 7 and nuclear recoil events can be achieved when the photo- 
electron statistics is larger than 20, which corresponds to an electron-equivalence energy 
threshold of about 5 keV, or 40—50 keV recoil kinetic energy, in the adopted detector 
configuration of 0.56 kg target mass. In realistic Dark Matter experiments, the modular 
mass for the CsI(Tl) target will have to be bigger, such as at the range of several kg. To 
maintain or even improve on such threshold, the light transmission within the crystal and 
the optical coupling between the crystal surface and the PMT photo-cathode will have 
to be optimized. Larger PMT readout surfaces as well as green-extended photo-cathodes 
to match the spectral emissions of CsI(Tl) can be used. Photon detectors with higher 
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quantum efficiencies, such as avalanche photo-diodes, can also be considered, though these 
devices tend to be limited in the active surface area. 

Although this article focuses on data with CsI(Tl) crystal scintillators for Dark Matter 
searches, the software techniques are readily applicable to other detector systems for 
other experiments where pulse shape analysis of individual events can provide useful 
information. A universal reference pulse profile was adopted in the present studies, one 
for recoil and one for 7 events, for the entire energy range of interest. This is justified in 
view of Figures EK&b such that the results in Figures EK&b are valid for CsI(Tl) crystals. 
In other detector systems where there may be a stronger energy dependence of the pulse 
shapes (as in the case for 7-events with Nal(Tl) crystals [H]), the PSD procedures can be 
further refined by having a different reference profile for every energy bin. The profiles 
can be derived from measurements with neutron beam and radioactive 7-sources like in 
this work. In addition, the conclusions on the relative merits among the different PSD 
methods are expected to be applicable to other pulse shape analysis problems where the 
statistics are marginal. 

Besides differentiating /5/7-background from nuclear recoil events, these studies may 
help to lower the detection threshold by suppressing electronic noise and microphonism 
where the pulse shapes are in general different from those of the signals. Experiments 
which need both low threshold and background may potentially benefit from these tech- 
niques. Alongside with Dark Matter experiments, such requirements are critical in the 
search of neutrino magnetic moments jl] and in the measurement of the coherent scatter- 
ings of the neutrinos on the nuclei |2nj . 
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Crystal 


Event Type 


Rise Time 

[to (ns)] 


Decay Time Constant 


Ratio (r) 


Fast Comp. 


Slow Comp. 

[t2 (/is)] 


CsI(Tl) 


nuclear recoils 


280±50 


0.54±0.1 


2.0±0.2 


0.29±0.02 


CsI(Tl) 


7 


296±50 


1.3±0.1 


4.50±0.4 


0.58±0.06 


Csl(pure) 


7 


~0.55 


0.19±0.02 







Table 1: Fitted rise and decay time constants as well as the ratio between slow and fast 
decay components for recoil and 7 events measured by CsI(Tl) and undoped Csl. 
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Figure 1: The average pulse shapes at 4.8 keV electron-equivalence energy for the nuclear 
recoil and 7-background events directly from the neutron beam measurements. 
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Figure 2: Comparison of average pulse shapes between 7- and nuclear recoil events: (a) 
recoil pulse at 4.8 keV electron-equivalence energy as compared to 7-events from 5 keV 
to 40 keV; (b) 7-pulse at 20 keV as compared to recoil profiles from 4 keV to 11 keV 
electron-equivalence energy. ^ ^ 
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Figure 3: Typical single (a) measured 7- and (b) simulated nuclear recoil events at 



Npe ~ 20. 
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Figure 4: Comparison of the distributions of (a) F(X) and (b) LR parameters between 
measured and simulated data denoted by solid circles and histograms, respectively. 
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Figure 5: The partial charge versus total charge at the high (MeV) energy range in 
a Csl(Tl) crystal, showing excellent (>99%) pulse shape discrimination capabilities to 
differentiate events due to a's and 7's. The a-events are from an ^"^^Am source (kinetic 
energy 5.49 MeV) placed on the surface of the crystal, while the 7-events are due to 
ambient radioactivity. 



19 




Figure 6: Typical separations of the (a) (t), (b) F(X) , and (c) LR parameters at Npe = 20 
between nuclear recoil (Dnr, in solid histograms) and 7 (D^, in dotted histograms) events 
with the mean time, neural network and likelihood ratio methods, respectively. 
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Figure 7: The variations of the figures of merit (a) ego and (b) /90 with Npe with the three 
different techniques apphed to simulated nuclear recoil and 7 data, respectively. Dotted 
lines indicate survival probabilities of 7 and recoil events in (a) and (b), respectively. The 
statistical uncertainties are smaller than the data points. 
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